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ABSTRACT 

We investigate the impact finite simulation box size has on the structural and 
kinematic properties of Cold Dark Matter haloes forming in cosmological simulations. 
Our approach involves generating a single realisation of the initial power spectrum of 
density perturbations and studying how truncation of this power spectrum on scales 
larger than L cut = 2n/k C ut affects the structure of dark matter haloes at 2 = 0. 
In particular, we have examined the cases of L cut = / cu tTbox with / cut = 1 (i.e. 
no truncation), 1/2,1/3 and 1/4. In common with previous studies, we find that the 
suppression of long wavelength perturbations reduces the strength of clustering, as 
measured by a suppression of the 2-point correlation function £(r), and reduces the 
numbers of the most massive haloes, as reflected in the depletion of the high mass 
end of the mass function n(M). Interestingly, we find that truncation has little impact 
on the internal properties of haloes. The masses of high mass haloes decrease in a 
systematic manner as L cut is reduced, but the distribution of concentrations is unaf¬ 
fected. On the other hand, the median spin parameter is ~ 50% lower in runs with 
/cut < 1- We argue that this is an imprint of the linear growth phase of the halo’s 
angular momentum by tidal torquing, and that the absence of any measurable trend 
in concentration and the weak trend observed in halo shape reflect the importance of 
virialisation and complex mass accretion histories for these quantities. These results 
are of interest for studies that require high mass resolution and statistical samples 
of simulated haloes, such as simulations of the population of first stars. Our analysis 
shows that large-scale tidal fields have relatively little effect on the internal properties 
of Cold Dark Matter haloes and hence may be ignored in such studies. 
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1 INTRODUCTION 

Over the last twenty years, the cosmological N-Body simu¬ 
lation has come to be firmly established as the preeminent 
tool for theoretical studies of th e formati on and evolution of 
structure in the Universe (see iBertschingeil 1 1998{l for a rel¬ 
atively recent review). Whereas early studies explored the 
gross features of different cosmological models and led to 
ergence of th e Co ld Dark Matter (CDM) paradigm 

et al. JJIIS; [ White, Frenk fe Dav^ lljai ; Davis et al] 
the most recent generation of simulations are be¬ 
ing used to provide detailed predictions of the CDM model 
and in particular its ACDM variant. These include the be¬ 
haviour of CDM halo mass profiles on scales that can be di- 
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rectly compared with observations (e.g . iHavas hi et al • 2003 : 
INa.varro et a.1.1 120041 : iDiema.nd et a,l. I 120041) and sonhisti- 
cated semi-analytical models of galaxy formation built upon 
merger trees harvested from high resoluti on simu lati ons 
(c.f. t he Millennium Simulation described in ISnringel e/all 
120051) . The N-Body approach is not without its limitations, 
however, and given its significance, it is of fundamental im¬ 
portance to understand what these limitations are and how 
they affect the interpretation of the results of simulations. 

Such considerations generally lead to “convergence” 
studies that seek to quantify the degree to which the results 
of a simulation are affected by choices made in its sett in g up 
and r unnin g. For examp l e, bot h INa.varro. Frenk fc W hite 
lll99fi!) and iMoore et all (1998) demonstrated that their 
results on the form of CDM halo mass profiles were not 
compromised by forc e reso lu tion (or “soften in g”) a s well 
as starting redshift llNavarro. Frenk fc White Ill996 ) and 
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pa rticle numbe r ( [Moore et alJ 1199 8.1 . The recent study 
of IPower et al • l2003h established a set of convergence 
criteria (reaffirmed hv IHa.vashi et all 120041 1 that determine 
the radial range over which the circular velocity profile 
of a simulated CDM halo is reliable to better than 10% 
accuracy, for a number of simulation parameters including 
softenin g, time integration accuracy, and particle number 
(see also|Diemand^et_al . 20041 ; lBinnevl2004ll . More generally, 
iKnebe et alJ (feOQojf investigated how properties such as the 
correlation function, mass function and central densities 
of dark matter haloes varied with softening and time 
integration accuracy at fixed particle mass in simulations 
run with different N-Body codes. 

Recently, iBagla fc Rav fcoOS . hereafter BR05) have 
demonstrated that size of simulation box Lbox can also af¬ 
fect certain properties of the dark matter halo population, 
in particular the 2-point correlation function £(r) and the 
mass function n(> M). Lbox sets the longest wavelength 
perturbation that can be resolved in a simulation and it 
follows that reducing Lb ox suppresses the contribution of 
large scale perturbations to the power spectrum. The effect 
of this suppression was highlighted bv lGelb fc Bertschingeii 
Jl994bl. hereafter GB94), who examined the importance of 
long wavelength perturbations in the initial conditions for 
£(r) and the correlation length ro in the context of the 
SCDM model. They investigated how ag (the linear mass 
variance in a sphere of radius 8/i _1 Mpc) and consequently 
£(r) varied with Lbox (see their figures 1 and 2), and found 
that Lbox > lOOMpc (or 50h _1 Mpc for h = 0.5) was re¬ 
quired to correctly recover both ag and the amplitude of 
£(r). In other words, studies that wish to accurately charac¬ 
terise the clustering properties of dark matter haloes require 
large simulation boxes. 

How large the simulation box must be to accurately re¬ 
cover £(r) was the question addressed by BR05. They noted 
that finite box size can affect not only £(r) but also the high 
mass end of the mass function n(> M); as Lb ox is reduced 
below some threshold value, the numbers of the most mas¬ 
sive haloes decrease in a systematic fashion . The m ass func¬ 
tion n (> M ) can be characterised by the lPress fc Schechteil 
(1974:) form which is a relatively simple function of the lin¬ 
ear mass variance a(M). For a given mass M = 4n/3 pR 3 
with p mean density, 

/> OO 

a 2 (M) = / P(k)W 2 (kR)d 3 k, (1) 

J2ir/L box 

where P(k) is the linear power spectrum, W 2 (kR) is the 
top-hat filter, and the lower limit of the integral 2n/Lbox 
corresponds to the fundamental mode in the simulation 
box. BR05 investigated how varying Lb ox impacted on the 
form of n(> M) and used this as a simple criterion for 
determining how large Lbox must be to recover n(> M) 
(and consequently £(r)) to a given accuracy. 

The studies of GB94 and BR05 clearly indicate that 
large simulation boxes are required if we wish to recover 
accurate mass and two point correlation functions of the 
dark matter halo population, but does the choice of Lbox also 
affect the internal properties of the haloes? This question is 
of interest for studies that require the high spatial resolution 
offered by simulations of small volumes but that may not 


require accurate clustering information, suc h as stu dies of 
the first objects at high redshift (e.g. lAbel et al . 20(ill . and 
it is one that we shall address in the present study. 

Why might we expect finite Lb ox to be of importance 
for the internal structure of haloes? The suppression of 
long wavelength perturbations will modify the “global” 
tidal field that a dark matter halo is subject to, which 
may have an impact on the halo’s shape, its angular 
momentum content (e.gJWiite|^)84) and the infall patter n 
of substructures (e.g. T/GiberteTaLl 120041 : iBensonl l2005lf . 
Similarly, the linear mass v ariance controls the formation 
time of dark matte r haloes ([Ncwarro, F renk fc Whit,ellT99fil : 

iBullock et al l200ltlEke et alfeoQl ) and so we might expect 

haloes to form at systematically later times in smaller 
boxes with lower “effective” values of ag; this may then 
affect both the mass of the halo and its central density (or 
concentration) measured at 2 = 0. We also note that small- 
and large-scale modes are coupled during the non-linear 
clustering phase and there is a tran sfer of power from large 
scale s down to small scales (e.g. IBagla fc Padmanabhanl 
Il997f) : neglecting this power may leave an imprint on the 
formation and evolution of gravitationally bound objects. 
On the other hand, dark matter haloes are, by definition, 
virialised systems and much of the information that was 
present during the initial stages of their collapse will be 
er ased during vi rial isati on (c.f. t he un iversal mass profile 
of [Navarro. Frenk fc Whitel 11996) . Il997l) , so it is not clear 
whether a finite Lb ox will have an effect at all. Nevertheless, 
it is important to investigate this question and determine 
how severe a limitation it might be. 

The outline of this paper is as follows; in the next sec¬ 
tion, we briefly describe the simulations we have used in 
this study before presenting our results (§ m We confirm 
the findings of GB94 and BR05 (§ 13.211 before examining 
how various characteristic properties of the halo population 
- concentrations, shapes and spins - are affected by the sup¬ 
pression of long wavelength perturbations on scales larger 
than L cut = f cut Lbox- Finally, we discuss our results in §0 
and offering some concluding remarks in §[£1 


2 THE SIMULATIONS 

We have run a set of cosmological N-body simulations fol¬ 
lowing the formation and evolution of structure in a set of 
simulation boxes of side Lbox = 128/i _1 Mpc, assuming a flat 
cosmology with flo = 0.333, 12a = 0.667, h = 0.667 and ag = 
0.88. Each simulation is single-mass and contains 256 3 par¬ 
ticles, implying a particle mass of m v = 1.15 x IO^/i^Mq; 
we adopt a starting redshift of 2 a tart = 40.0 in each case. 
For each of the runs, we use the same realisation (i.e. the 
same amplitudes and phases) of the power spectrum 1 but 
we truncate the power on scales k ^ fc cu t = 2n/L cu t, where 
L cut is a fraction / cut of Lb ox and / cu t = 1 (i.e. no trun¬ 
cation), 1/2, 1/3 and 1/4. These correspond to cut-offs of 
L cu t = 128,64,43 and 32 h -1 Mpc or fc cut ~ 0.05,0.1,0.15 
and 0.2 h Mpc -1 respectively. Truncating the power spec¬ 
trum in this way mimics the effect of varying Lbox while 

1 Generated using CMBFAST ISeliak fc Zaldarriaga| [l99(Jf 
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allowing the impact of the truncation on individual haloes 
to be investigated. 

The simulations were performe d usi ng the parallel 
TreePM code GADGET2 llSpringel et alJEoOSl) with a fixed co¬ 
moving softening of 0.01h^ 1 Mpc and a timestep parameter 
of Tj = 0.02. 

Groups were identified using AHF 2 , a mod ifica tion of 
the MHF algorithm that was presented in iGill et al.l <2004 . 
Haloes are located as peaks in an adaptively smoothed den¬ 
sity field of the simulation using MLAPM’s grid hierarchy and 
a refinement criterion that matches the force resolution of 
the actual simulation carried out with GADGET2 (i.e. 5 par¬ 
ticles per cell); local potential minima are computed for 
each of these peaks and the set of particles that are grav¬ 
itationally bound to the halo are returned. For every halo 
(either host or satellite) we calculate a suite of canonical 
properties (e.g. velocity, mass, spin, shape, concentration, 
etc.) based upon the particles within the virial radius. The 
virial radius R v n is defined as the point where the den¬ 
sity profile (measured in terms of the cosmological back¬ 
ground density pb) drops below the virial overdensity A v i r , 
i.e. M(< J? v i r )/(47rf?5 ir /3) = A v ; T pb- This threshold A v n 
is based upon the dissipationless spherical top-hat collapse 
model and is a function of both cosmological model and 
time. For 2 = 0 it amounts to A v n = 340. 


3 RESULTS 

We now present the results of our analysis. Firstly, we give 
a visual impression of the impact of truncating the power 
spectrum at large scales on the overall clustering pattern 
within the simulation volume. Secondly, we consider some 
bulk measures of the dark matter halo population, including 
the mass function n(> M), 2-point correlation function £(r) 
and power spectrum P(k). Finally we examine the internal 
properties of haloes, including concentrations, shapes and 
spin parameters. 

I 11 what follows, we compare and contrast the runs in 
which the initial power spectrum is truncated with the fidu¬ 
cial run ChopNo, whose initial power spectrum was not trun¬ 
cated (i.e. /cut = 1). The “truncated” runs are called Chop64 
(i.e. /cut = 1/2), Chop43 (i.e. / cu t = 1/3) and Chop32 (i.e. 
/cut = 1/4), respectively. 

3.1 Spatial Distribution at 2 = 0 : A Visual 
Impression 

Figure Q shows projected maps of the dark matter density 
distribution within a thin slice taken through each of the 
four simulations at 2 = 0; for clarity, we have chosen a slice 
thickness of 10% the box length (i.e. 12.8/i _1 Mpc) taken 
through the centre of the box. 

Although it is possible to pick out features that are 
common to all four runs, the magnitudes of the correspond¬ 
ing overdensities decrease as L cu t is decreased. The network 
of filaments in the fiducial run is more striking than in the 
“truncated” runs - that is, the fraction of volume that is 

2 AHF can be downloaded from the following web page 
http: //ww. aip. de/People/aknebe/AMIGA 
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Figure 2. The Cumulative Mass Function n(> M). We show 
the number of haloes with mass greater than M in each of the runs 
at redshifts z = 5, 3,1 and 0. The heavy solid curves correspond 
to the fiducial run while the dotted, dashed and dotted-dashed 
curves correspond to the runs set up with truncated power spec¬ 
tra. We also plot the predicted Press-Schechter mass function for 
the appropriate redshift (light solid curves). See text for further 
details. 

at or above a given overdensity threshold is smaller in the 
runs in which large scale power was suppressed The central 
overdensities of individual haloes are comparable in the re¬ 
spective runs, but these haloes tend to occupy regions that 
are less overdense in the “truncated” runs. Similarly, the 
local mean number density of haloes appears lower in the 
“truncated” runs. Both of these latter observations are con¬ 
sistent with the GB94 and BR05 results that finite box size 
impacts on £(r) and consequently the correlation length ro. 

3.2 Bulk Properties of the Halo Population 

The projected maps of the dark matter density distribution 
suggest that haloes are less clustered and that there are 
fewer massive haloes in the “truncated” runs. We now 
quantify these differences by studying the halo mass 
function n(> M), the correlation function £(r) and the 
non-linear power spectrum P{k). 

We begin by comparing cumulative halo mass functions 
n(> M) for each of the runs in figure [ 5 ] these mass func¬ 
tions have been constructed from the AHF halo catalogues 
at redshifts 2 = 0,1, 3 and 5. The heavy solid curves corre¬ 
spond to the fiducial ChopNo run while the dotted, dashed 
and dotted-dashed curves correspond to the Chop64, Chop43 
and Chop32 runs respectively. For reference, we also plot 
the appropriate Press-Schechter mass function for the given 
redshift (light solid curve in each of the panels). 

This figure demonstrates that truncating the power on 
large spatial scales leads to a suppression of the number of 
high mass systems that form, and the effect is most notice¬ 
able in the runs where the truncation is most severe. We 
note that the effect is most pronounced at 2 = 0 but it is 
apparent at all redshifts. 

We investigate this effect further in figure [ 3 ] where we 
cross correlate the masses of individual haloes that have 
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Chop64 

Chop43 
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Figure 1. The projected density distribution in 12.8h _1 Mpc slices taken through the centres of each of the boxes. We have smoothed 
the particle mass using an adaptive Gaussian kernel and projected onto a mesh. Each mesh point is weighted according to the logarithm 
of its projected surface density, and so the “brighter” the mesh point, the higher the projected surface density. 


been matched in the fiducial and truncated runs, i.e. we 
plot the mass ratio of corresponding haloes against the mass 
in the fiducial run. This is achieved by identifying its 10% 
most bound particles in the ChopNo model and locating those 
haloes in the truncated runs that contain the largest fraction 
of these particles; this allows for unique matching. 

Here we observe a clear trend - the masses of the high 
mass systems are systematically smaller in the truncated 
runs relative to the fiducial run, and the effect becomes 
more pronounced as the truncation is applied more ag¬ 
gressively (compare the upper and lower panels). However, 
as we go to lower masses (10 12 — 10 13 /j -1 Mq), we find 
that the median halo mass is unaffected, although we 
note that the mass of any given halo can change by as 


much as a factor of 10 (with a rms variation of a factor of 2). 

Figure Q] shows how clustering as measured by the 
correlation function £(r) varies between different runs. The 
heavy lines correspond to the 850 most massive haloes 
within the simulation volume (according to a number 
density of roughly 4 x 10 -4 /i 3 Mpc -3 ), while the light 
lines correspond to the dark matter. Solid, dotted, dashed, 
dotted-dashed curves represent the fiducial runs and the 
truncated runs respectively. The upper panel shows £(r) 
while the lower panel shows the ratio of £(r) measured for 
each of the truncated runs relative to the fiducial run. The 
first point to note is that the haloes are more clustered than 
the dark matter - a well-known result. The second point is 
that the amplitude of £(r) for both the dark matter and 
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Figure 3. Comparison of Halo Masses in the truncated and fidu¬ 
cial runs. See text for further details. 


the haloes decreases as the cut-off scale L cu t decreases; the 
extent of the effect is clearly illustrated in the lower panel. 

As a final measure of the impact of truncating the ini¬ 
tial power spectrum, we measure the non-linear dark matter 
power spectrum at z = 0 for each of the runs, shown in fig- 
ure|3 At small scales there is good agreement but deviations 
from the fiducial run become apparent at k ~ lO/iMpc^ 1 , 
at which point the amplitude of the power spectrum is sup¬ 
pressed in the truncated runs relative to the fiducial run. 
This is clear in the lower panel where we show the ratio of 
the power spectra to the fiducial case. 


3.3 The Internal Properties of the Halo 
Population 

We have so far considered bulk properties of the halo 
population and shown that our results are consistent with 
those of GB94 and BR05. However, our primary motivation 
for this study is to determine whether the choice of box 



0.01 0.10 1.00 10.00 
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Figure 4. The 2-Point Correlation Function £(r) (upper panel) 
and the ratio of Cchopxx(r) in the truncated runs relative to the 
fiducial run £chopNo ( r ) (lower panel) for the dark matter (light 
curves) and the haloes (heavy curves). 


size influences the structural and kinematic properties of 
haloes, and this is the subject of this section. To this extent 
we are focusing on three relevant integral properties of dark 
matter haloes, namely concentration, spin parameter, and 
triaxiality. 


We begin by investigating the impact of the finite box 
size on concentration, which provides a measure of the cen¬ 
tral density of a halo. Typically a concentration c v ir is de¬ 
fined; 

Cvir = — (2) 

r s 

where r s is a “scale” radius (c.f. lNa.va.rro. Frenk fc Whitel 
119971 defined as the radius at which the logarithmic slope 
of the de nsity profile is —2) . In this study, we have adopted 
Ci /5 (e.g. lColin et alJfeoQCl iKnebe et alll2002l) . defined as 


„ R{l/5M viI ) 
°l/5 - -5- 

itvir 


(3) 


R(l/5 M v i r ) is the radius enclosing 1/5 the virial mass of the 
halo. This is a convenient, robust and model-independent 































6 Power & Knebe 



0.1 1.0 10.0 100.0 
k \h Mpc -1 ] 


Figure 5. The Dark Matter Power Spectrum P(k) (upper panel) 
and the ratio of power in a given mode in the truncated runs 
relative to the fiducial run (lower panel). 

measure of concentration, although it is a non-trivial func¬ 
tion of c v ir; for example, c v n = 5 corresponds to C 1/5 = 5 
while Cvir = 10 is approximately C 1/5 = 7. However, this is 
unimportant because we are interested in relative variations 
in concentration. 

We argued in § Q that the suppression of long 
wavelength perturbations affects the linear mass variance 
a(M), which should affect when the halo forms and hence 
what its concentration at z = 0 will be. Furthermore, 
previous st udies have shown that c v u i s a function of 
halo ma ss llNavarro. Fr enk fc WhiteHlifflit lEke et all 1200 it 
iBullock et al.ll20(nl) ' an d the density of the environment in 
which the halo resides iBullock et al.ll200ll hereafter B01). 
It follows that Cvir might be influenced by the finite box size; 

• The relation between halo mass and c v ir is such that 
Cvir is a decreasing function of increasing halo mass. Higher 
mass systems tend to form at later times when the mean 
density of the Universe was lower, and so the central den¬ 
sities and hence concentrations of these haloes are smaller. 
Suppressing the contribution from long wavelength pertur¬ 
bations results in a lower effective as, which should delay 
the formation time of haloes and hence reduce c v ir- We ex- 



Figure 6 . Distribution of Concentrations C 1 / 5 . Solid, dotted, 
dashed and dotted-dashed histograms correspond to the ChopNo 
(fiducial), Chop64, Chop43 and Chop32 runs respectively. The thin 
solid lines are derived by fitting the distribution to a Gaussian. 

pect this effect to be more dramatic for higher mass systems 
(c.f. figure [U. 

• Haloes in less dense environments tend to have lower 
concentrations than haloes of similar mass in higher density 
environments (B01). Figure^Jshows that both the dark mat¬ 
ter and haloes are less clustered in the truncated runs and 
thus reside in lower density environments relative to their 
counterparts in the fiducial run. 

Although there are two competing effects, we note that 
B01 observed that the trend with environment is stronger 
than the trend with mass. Moreover, it is the masses of the 
higher mass systems (> 10 13 /i -1 M@) that are most affected, 
where the effect on clustering affects both the dark matter 
and the haloes in a similar manner. Based on these observa¬ 
tions, we might expect a shift towards higher concentrations. 

In figure^ we show the distributions of C 1/5 for the 
four runs. Each peaks at C 1/5 ~ 5 (c v ir — 5), and while there 
are small differences in amplitude and offset in the position 
of the peak, there is no indication that the differences be¬ 
tween the distributions are significant. This is confirmed by 
computing the KS statistic for the truncated runs compared 
with the fiducial run; the probability that the differences 
are significant is small (<) 10 %) and negligible compared to 
the variation between runs with different realisations of the 
(untruncated) power spectrum. Restricting our analysis to 
only those haloes more massive than 10 13 /i _1 Mq does not 
affect the result. 

We checked for differences in both the mass accretion 
histories and the mean local density of haloes (c.f. figure 6 
of B01) between the truncated and fiducial runs, but again 
the differences are negligible. Haloes in each of the runs 
accrete similar fractions of their final mass between z = 1 
and z = 0 (e.g. ~ 50% for 10 12 h _1 Mq haloes) and the mean 
local density with a sphere of lh -1 Mpc centred on each 
halo in the different runs correlates with C 1/5 as we would 
expect. It seems that the concentration of individual haloes 
are sensitive to the suppression of long wavelength pertur¬ 
bations in the initial conditions (as we shall see in § 1 . 1 , 4ll . 
but those of a “typical” halo do not. This result is further 
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Figure 7. Comparison of halo concentrations C 1/5 in the trun¬ 
cated and fiducial runs. 
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Figure 8. Distribution of Spin Parameters in the fiducial (solid 
curve and histogram) and truncated (dotted, dashed and dotted- 
dashed curves and histograms) run s. Curves represent fits to the 
lognormal distribution favoured bv lBullock et al . I 2001 al) . 


The shapes and amplitudes of the distributions are 
all well described by the BOla lognormal functional form. 
And as the truncated models only display a marginal shift 
towards smaller spin parameters in figure [S] the effects 
of the (missing) large-scale modes are better revealed in 
figure 0 which can be considered an analogue to figure [3 
again; we match haloes in the truncated runs with their 
counterparts in the fiducial run that share the largest 
fraction of the 10% most bound particles of the halo. We 
can now see that there is an offset of ~ 50% between the 
fiducial run and the truncated runs, independent of halo 
mass; that is, the median spin is ~ 50% larger in the 
fiducial run for haloes of all masses. We also note that the 
spin of an individual halo can vary by as much as a factor 
of 10 (rms variation of a factor of ~ 3) between truncated 
and fiducial runs. 


confirmed when plotting the analogue to figure [3 i.e. the 
cross-correlation between the individual halo concentrations 
as shown in figure Q We again observe a scatter about 
the mean ratio, but on average this ratio is found to be unity. 

We have concentrated on the mass profile so far, but 
we also argued that the angular momenta and the shapes of 
haloes should be affected by the suppression of long wave¬ 
length perturbations in the initial conditions. Figure|H]shows 
the distribution of spin parameters A' in the truncated and 
fiduci al runs. Note that we have adopted the lBullock et all 
fcOOlal. hereafter BOla) definition in this figure; 

A' = —=^~ -. (4) 

V2 MVR 

The light curves correspond to fits to the lognormal distri¬ 
bution advocated by BOla 3 . 

3 We confirm that these distributions are in good agreement with 
those derived for spins defined by A = J\E\ 1 / 2 /GM 5 / 2 . 


Finally, we examine the shapes of haloes forming in the 
truncated runs and compare with those that form in the 
fiducial run. For convenience, w e compute the triax iality pa¬ 
rameter that was introduced bv lFranx et al. . f 1991 1. defined 
as 



(5) 


where a ^ b ^ c are the lengths of the major, intermedi¬ 
ate and minor axes, respectively, as defined by all particles 
interior to R v i T . The resulting distributions are shown in fig¬ 
ure E3 We observe again a trend for a marginal decrease 
in triaxiality when reducing the power on large scales. This 
is further confirmed by figure fTTI where we plot the ratio of 
triaxiality for all cross-identified halos. The signal is similar 
(yet weaker) to the one already reported for the spin param¬ 
eter: the more large-scale power has been ignored the smaller 
the triaxiality. However, it is not as pronounced, especially 
not for the ratio of ChopNo to Chop32. 
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Table 1. Physical Properties of Halo. We show the virial mass, M v j r ; the virial radius, -R v ir! the radius at which the circular velocity 
peaks, i?max! the peak circular velocity, W,max5 the NFW concentration, c v ir; the spin parameter, A; and the axis ratios, b/a and c/a. 


/cut 

M vir 

[lO^fc-iMo] 

R-vir 

[/i _1 Mpc] 

Rmax 

[/i -1 Mpc] 

V c,max 

[kms -1 ] 

Cvir 

A 

( b/a,c/a ) 

T 

1 

5.336 

1.65 

0.32 

1384 

8.3 

0.0089 

(0.73,0.43) 

0.57 

1/2 

2.625 

1.30 
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Figure 9. Comparison of Halo Spin Parameters in the fiducial 
and truncated runs. See text for further details. 


3.4 Case Study 

We have found that the internal properties of haloes are ro¬ 
bust in the sense that the concentrations and widths of their 
distributions and the values of their medians are practically 
unaffected as the “effective” box size is varied. However, we 
have also found that properties such as the triaxiality T and 
spin X' of individual haloes can vary by a factor of several 



Figure 10. Distribution of Triaxialities in the fiducial and trun¬ 
cated runs. See text for further details. 


between those truncated runs and the fiducial run. In this 
short case study, we illustrate how changing the effective 
box size impacts on an individual halo. 

We considered a number of characteristic internal prop¬ 
erties of a single well resolved halo that we have identified 
in each of the runs; these are summarised in table Q while 
projections of the particle distribution within 1.5 times the 
virial radius are shown in figure [H] We have chosen the most 
massive halo in the fiducial simulation for our analysis be¬ 
cause we expect the differences to be most pronounced on 
this mass scale. 

Inspection of table Q reveals that truncation of the 
power spectrum has a dramatic effect on the halo’s mass, 
as expected from figure [3 - a factor of ~ 2 between the 
fiducial run and the three truncated runs. This is reflected 
in the decrease in concentration (by 30 — 40%). Interest¬ 
ingly, the spin parameter increases steadily by up to a fac¬ 
tor of ~ 2 between the fiducial and the Chop32 truncated 
run. Also, we note that the semi-minor axis c/a remains 
fixed (at ~ 0.43) but the semi-major axis steadily decreases 
(from 0.73 to 0.46), leading to an increase in triaxiality (from 
T=0.57 to T=0.98). 


4 DISCUSSION 

The results presented in the previous show that the size of 
simulation box has little impact, if any, on the internal struc¬ 
ture and kinematics of simulated Cold Dark Matter haloes. 
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Figure 11. The ratio of the triaxiality parameter T in cross- 
identified halos. 


At first, this result is a little surprising, for reasons given in 
§ 0 Although we focused on concentration, shape and spin 
parameter because we expect these quantities to be most 
susceptible to the contribution of long wavelength pertur¬ 
bations, we also examined a number of other quantities, in¬ 
cluding radial velocities, velocity anisotropies and the virial 
ratio 2T/|VF|. In all cases, the distributions recovered from 
the truncated runs are indistinguishable from those derived 
from the fiducial run. 

However, it can be argued that this result is not so sur¬ 
prising. The impact of large scale perturbations on the value 
of <T8 noted by GB94 is certainly noticeable - of order a 40% 
reduction for a box size of Abox = 32A -1 Mpc relative to the 
“true” value - but the effect is less pronounced as we probe 
smaller mass scales. For example, the difference is less than 
~ 10% at 2/i 1 Mpc for the same box size. In other words, 
large scale perturbations become less important as we probe 
smaller masses. We can understand this in figure [3 which 
both confirms the result of BR05 that it is the higher mass 
systems that are most affected by the truncation of long 



-5 0 5 -5 0 5 

[Mpc/h] 


Figure 12. Projected particle distribution in the neighbourhood 
of a 10 14 /i —1 Mq cluster mass halo in each of the four runs. The 
virial mass of the halo decreases as the power spectrum is trun¬ 
cated at progressively shorter wavelengths. 


wavelength perturbations and also highlights that the typi¬ 
cal or median halo mass below 10 13 /) -1 Mq is unaffected 
by truncation. 

Perhaps more pertinent ly, the ex istenc e of a universal 
mass profile llNavarro. Frenkfc Whitdtl99A (T9971 that pro¬ 
vides an adequate characterisation of haloes in a range of 
hierarchical cosmologies implies that large scale perturba¬ 
tions have little effect on the inner structure of dark mat¬ 
ter haloes. Indeed, truncating the power spectrum at short 
wavelengths prod uces a ha lo mass pro file that is w ell de¬ 
scribed bv the iNava.rro. Frenk White! dl 99(11199711 form 
iMoore et allMiffilKnebtH^rLl^OO^h Nevertheless. our re¬ 
sults show that the internal structure of haloes (in the spher¬ 
ically averaged sense) are unaffected by the introduction of 
a long wav elength cut-off in the p ower spectrum (see also 

Similarly, the weakness of the trend in the distribution 
of halo shapes most likely reflects the relevance of the de¬ 
pendence of shape on merging history. While the infall of 
material on the halo is sensitive to the global tidal field, 
and therefore long wavelength perturbations, a change in 
the tidal field does not necessarily suppress infall; rather it 
will change the time at which it occurs and therefore will 
not affect halo shape in an average sense. 

The small effect noted in the spin distributions - an off¬ 
set in the median values of the truncated and fiducial values 
- suggests that long wavelength perturbations play some role 
in shaping the spin and consequently the angular momen¬ 
tum content of the halo. Tidal torques make an important 
contribution to the acquisition of angular momentu m of the 
halo during its growth in the linear regime ijWhitdil 984i l but 
we do not expect the halo’s subsequent non-linear evolution 
to depend on the large scale tidal field. The small differ¬ 
ence likely represents an imprint of the angular momentum 
acquired during the haloes linear growth. 
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5 CONCLUSIONS 

The aim of this study has been to establish whether finite 
box size has a measurable effect on the internal properties 
of simulated Cold Dark Matter haloes. The results of our 
analysis suggest that the effects, if present, are small and 
not statistically significant in most cases. Of the principal 
quantities we have examined - concentration, shape and the 
spin parameter - we find that spin shows the most prominent 
effect; the median spin parameter is ~ 50% smaller in our 
truncated runs, independent of the mass of the system. We 
argue that this is an imprint of the linear growth phase of 
the halo’s angular momentum by tidal torquing, and that 
the absence of any measurable trend in concentration or 
strong trend in shape reflect the importance of virialisation 
and complex mass accretion histories for these quantities 
respectively. 

These results are useful because they clarify what prop¬ 
erties of simulated dark matter haloes are affected by finite 
box size, and the severity of these effects. Indeed, they are of 
some importance because they reveal that studies of the in¬ 
ternal properties of statistical samples of haloes that do not 
require clustering information that are based on high resolu¬ 
tion simulations of small volumes are not compromised. The 
abundance of high mass systems will be underestimated, 
but the properties of low- to intermediate-mass systems are 
reliable. This is of interest to studies of, for example, the 
formation of the fi rst generat ion of population III stars at 
high redshift fe.g. lAhel et a.lil2002 ); in these cases, cluster¬ 
ing information is not as important (although it would be 
for reionisation studies, say). 

In conclusion, our study has demonstrated that the ef¬ 
fects of finite box size are negligible for the internal proper¬ 
ties of dark matter haloes. 
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